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1 Introduction 

Consider a manifold constructed by identifying the boundaries of Euclidean tri- 
angles or Euclidean tetrahedra. When these form a closed topological manifold, 
we call such spaces piecewise fiat manifolds (see Definition [l| as in [5] . Such 
spaces may be considered discrete analogues of Riemannian manifolds, in that 
their geometry can be described locally by a finite number of parameters, and 
the study of curvature on such spaces goes back at least to Regge [35] . In this 
paper, we give a definition of conformal variation of piecewise flat manifolds in 
order to study the curvature of such spaces. 

Conformal variations of Riemannian manifolds have been well studied. While 
the most general variation formulas for curvature quantities is often complicated, 
the same formulas under conformal variations often take a simpler form. For 
this reason, it has even been suggested that an approach to finding Einstein 
manifolds would be to first optimize within a conformal class, finding a mini- 
mum of the Einstein-Hilbert functional within that conformal class, and then 
maximize across conformal classes to find a critical point of the functional in 
general [T\. Finding critical points of the Einstein-Hilbert functional within a 
conformal class is a well-studied problem dating back to Yamabe j48j , and the 
proof that there exists a constant scalar curvature metric in every conformal 
class was completed by Trudinger [45 , Aubin [2\ , and Schoen flD] (see also [2^ 
for an overview of the Yamabe problem). 

Implicitly, there has been much work on conformal parametrization of two- 
dimensional piecewise fiat manifolds, many of which start with a circle packing 
on a region in or a generalized circle packing on a manifold. Thurston 
found a variational proof of Andreev's theorem ([ll] and conjectured that 
the Riemann mapping theorem could be approximated by circle packing maps. 
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which was soon proven to be true [32] (see also for an overview of the 
theory). Another direction for conformal parametrization appears in |38j . |28j . 
and |42j . These and other works produce a rich theory of conformal geometries 
on surfaces and have led to many beautiful results about circle packings and 
their generalizations. The theory developed in this paper unifies several of 
these seemingly different notions of conformality to a more general notion. It 
also allows an explicit computation of variations of angles which allows one to 
glean geometric information. The geometric interpretation of the variations of 
angles was known in some instances (e.g., [23] [17]), but the proofs were explicit 
computations, which made them difficult to extend to more general cases. One 
of the main contributions of this paper is to show how these computations may 
be done in a more simple, geometric way which easily generalizes. 

Existing literature on conformal parametrization of three-dimensional piece- 
wise flat manifolds is much more sparse. A notion was given by Cooper and 
Rivin |12j which takes a sphere packing approach, and a rigidity result was pro- 
duced (see also [37| and [H]). However, this theory requires that edge lengths 
come from a sphere packing, which is a major restriction of the geometry even 
on a single tetrahedron. In [T7], the author was able to show by explicit com- 
putation that the variations of angles are related to certain areas and lengths 
of the piecewise flat manifold (in actuality, one needs the additional structure 
of a metric as described below) . The theory developed in this paper generalizes 
this result to a general class of three dimensional piecewise flat manifolds. This 
generalization allows a geometric understanding of the variation of angles in a 
three-dimensional piecewise flat manifold under conformal variations, and the 
space of conformal variations is quite large and need not depend on the initial 
distribution of the edge lengths (unlike |17j , where one must assume that the 
metric comes form a sphere packing structure) . 

The variation formulas for the curvature allow one to introduce a theory 
of functionals closely related to Riemannian functionals such as the Einstein- 
Hilbert functional. In two dimensions, many of these functionals are well stud- 
ied, originally dating back to the work of Colin de Verdierc [11 . In dimensions 
greater than three, the generalization of the Einstein- Hilbcrt functional was 
suggested by Regge [3S] and has been well studied both in the physics and 
mathematics communities (see j22j for an overview). Recently, the functional 
was used to provide a constructive proof of Alexandrov's theorem that a surface 
with positive curvature is the boundary of a polytope [1]. In this paper, we give 
a general construction for two-dimensional functionals arising from a conformal 
structure. We also consider variations of the Einstein-Hilbert-Regge functional 
with respect to conformal variations. Variation of this functional gives rise to 
notions of Ricci flat, Einstein, scalar zero, and constant scalar curvature met- 
rics on piecewise flat manifolds. Our structure allows one to consider second 
variations of these functionals around fixed points, and give rigidity conditions 
near a Ricci flat or scalar zero manifold. An eventual goal is to prove theorems 
about the space of piecewise flat manifolds analogous to ones on Riemannian 
manifolds, for instance [26] [3Tj . 

Certain curvatures considered here have been shown to converge in measure 
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to scalar curvature measure by Cheeger-Miiller-Schrader in [5]. The proof in 
the general case does not appear to give the best convergence rate, and it is an 
open problem what this best convergence rate may be. It would be desirable to 
have a more precise control of the convergence and to prove a convergence of 
Ricci curvatures or of Einstein manifolds on piecewise flat spaces to Riemannian 
Einstein manifolds. Although the convergence result shows convergence to scalar 
curvature measure, it has been suggested that these curvatures are analogous 
to the curvature operator on a Riemannian manifold [7]. 

This paper is organized as follows. Section [2] gives definitions of geomet- 
ric structures on piecewise flat manifolds in analogy to Riemannian manifolds 
and shows the main theorems on variations of curvature functionals. Section |3] 
derives formulas for conformal variations of angles. Section [4] translates these 
results to variations of curvatures and curvature functionals. Section [s] discusses 
some of the conformal structures already studied and shows how they fit into 
the framework developed here. Finally, Section [6] discusses discrete Laplacians, 
when they are negative semidefinite operators, and how this implies convexity 
results for curvature functionals and rigidity of certain metrics. The main the- 
orems in the paper are Theorems 29 and[3l]on the variations of angles, which 
could easily be applied to extend these results to the case of manifolds with 
boundary. Theorems |32| a nd |34| on the variation of curvature, which give ana- 
logues of the variation of scalar curvature under conformal deformation of a 
Riemannian metric. Theorems [23] and [24] on the variation of curvature function- 
als. Theorems [40] and ]4T] on convexity of curvature functionals, and Theorems 



43 and 44 on rigidity of zero scalar curvature and Ricci fiat manifolds. 



2 Geometric structures and curvature 
2.1 Metric structure 

We will consider certain analogues of Riemannian geometry. A Riemannian 
manifold (M", 5) is a smooth manifold M together with a symmetric, positive 
definite 2-tensor g. A piecewise fiat manifold is defined similarly to the defini- 
tions in [S]. 

Definition 1 A triangulated manifold (M, T) is a topological manifold M to- 
gether with a triangulation T of M. A (triangulated) piecewise fiat manifold 
[M, T, i) is a triangulated manifold (M, T) together with a function i on the 
edges of the triangulation such that each simplex can be embedded in Euclidean 
space as a (nondegenerate) Euclidean simplex with edge lengths determined by 
£. 

Nondegeneracy can be expressed by the fact that all simplices have positive 
volume. This condition can be realized as a function of the edge lengths using 
the Cayley-Menger determinant formula for volumes of Euclidean simplices. 

In this paper we will consider only closed, triangulated manifolds although 
the definitions could be extended to more general spaces. We will describe sim- 
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plices as . . . ,k} , where k are natural numbers. The length associated 
to an edge {i,j} will be denoted £ij, area associated to {i,j, k} will be denoted 
Aijk, and volume associated to {i,j,k,£} will be denoted by Vijke- Note that 
once lengths are assigned, area and volume can be computed using, for instance, 
the Cayley-Menger determinant formula. We will also use the notation jijk to 
denote the angle at vertex i in triangle {i, j, k} , and sometimes drop jk when 
it is clear which triangle we are considering. A dihedral angle at edge in 
{i,j, k,£} will be denoted Pij^ki and k£ will be dropped when it is clear which 
tetrahedron we are considering. In all of the following cases, the indices after 
the comma will be dropped when the context is clear. 

Definition 2 Let V (T) denote the vertices of T, let E (T) denote the edges of 
T, and let (T) denote the directed edges in T (there are two directed edges 
(i, and associated to each edge {i,j}). For any of these vector spaces X, 
let X* space of functions h : X ^ R. 

Note that, for instance, if d G -B+ (T) then d = j)eE+ d(i,j) (pij, where 
(pij is the standard basis of E+ (V)* . We will use dij (as in Definition |5| to 
denote either d{i,j) or the function d(i, j) (pij, and similarly with elements of 
V{T)* and E{T)* . 

Remark 3 We are implicitly assuming that the list of vertices determines the 
simplex uniquely. This is just to make the notation more transparent. We 
could also have indexed by simplices, such as i„2,A^3,"f^0(-^3, etc. This latter 
notation is much better if one wants to allow multiple simplices which share the 
same vertices. 

Remark 4 A piecewise flat manifold is a geometric manifold, in the sense that 
it can be given a distance function in much the same way that a Riemannian 
manifold is given a distance function, i.e., by minimizing over lengths of curves. 

The definition ensures that each simplex can be embedded isometrically in 
Euclidean space. The image of vertex i in Euclidean space will be denoted Vi, 
the image of edge will be denoted w^Wj, etc. 

Piecewise flat manifolds arc not exactly the analogue of a Riemannian man- 
ifold we will consider. 

Definition 5 Let [M, T) be a triangulated manifold. A piecewise flat pre- 
metric is an element d G E^ [T) * such that (M, T, I) is a piecewise fiat manifold 
for the assignment £ij = dij + dji for every edge . A piecewise flat pre- 

metric d is a metric if for every triangle {i,j, k} in T, 

dl + d% + dl = d% + + d% . (1) 

A piecewise flat, metrized manifold {M,T,d) is a triangulated manifold {M,T) 
with metric d. 

For future use, we define the space of piecewise flat metrics on (M, T) . 
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Definition 6 Define the space met {M, T) to be 

met {M, T) = {de E+ (T)* : (Af , T, d) is a piecewise flat, metrized manifold^ . 

As shown in [T^, condition (jlj ensures that every simplex has a geomet- 
ric center and a geometric dual which intersects the simplex orthogonally at 
the center. This dual is constructed from centers. Given a simplex embedded 
into space as {vi,V2, ■ ■ ■ ,w„} , we have a center point to the simplex given by 
Ci23---n- This point can be projected onto the (n — l)-dimensional simplices and 
successively projected onto all simplices, giving centers Cy...fc for all subsets of 
{1, 2, . . . , n} . The centers can be constructed inductively by starting with cen- 
ters of edges at a point which is a (signed) distance dij from vertex i and 
dji from vertex j. Then one considers orthogonal lines through the centers, and 
condition ([l]) ensures that in each triangle, there is a single point where these 
three lines intersect, giving a center for the triangle. The construction may be 
continued for all dimensions, as described in |19j . 

For simplicity, let's restrict to n < 4. We will denote the signed distance 
between C1234 and Cijk by hijk,(. and the signed distance between Cijk and Cij 
by hij^k- The sign is gotten by the following convention. If C1234 is on the same 
side of the plane defined by ViVjVk as the tetrahedron V1V2V3V4, then hijk/ is 
positive, otherwise it is negative (or zero if the point is on that plane). Similarly, 
if Cijk is on the same side of the line defined by ViVj as ViVjVk within that plane, 
then hij^k is positive. Since it is clear that hij = dij, we will not use the former. 
The side ViVj is divided into a segment containing Vi of length dij and a segment 
containing Vj of length dji such that iij — dij + dji. It is easy to deduce that 
hij^k and hijk.i can be computed by 

T ^ik ^ij COS ^i,jk 

n-ij,k = : 

smjijk 

and 

, _ hij^e — hij_k cos Pij^ke 

nijk.i — : — 7, • 

smfiijM 

See [12] or |3] for a proof. Importantly, these quantities work for negative values 
of the d's and /I's. We will also consider the dual area Aij^ke of the edge {i, j} 
in tetrahedron {i,j,k,£} , which is the signed area of the planar quadrilateral 
ci234CijkCijCije, where i,i,k,£ are distinct. The area is equal to 

AijM — 2 {hij,khijkx + hij^ghiji^k) ■ 

These definitions of centers within a simplex induce a definition of geometric 
duals on a triangulation (see [12] for details). In particular, we will need the 
lengths or areas of duals of edges, defined in two and three dimensions as follows. 

Definition 7 Let {AP,T,d) be a piecewise flat, metrized manifold of dimension 
2. Then edge is the boundary of two triangles, say {i,j,k} and £}. 

The dual length £*j is defined as 

tij — hij^k + hij^i. 
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Note that the two triangles can be embedded in the Euclidean plane together, 
and i*j is the signed distance between the centers of the two triangles. 

Definition 8 Let (AP,T,d) be a piecewise flat, metrized manifold of dimension 
3. Then the dual length i*^ (which is technically an area) is defined as 

k,e 

kj 

where the sum is over all tetrahedra containing the edge 

Notation 9 Most sums in this paper will be with respect to simplices, so a sum 
such as the one in Definition^ means the sum over all tetrahedra {i,j^k,f\ 
containing the edge , not the sum over all values of k and i (which would 

give twice the aforementioned sum). 

The dual length is the area of a (generalized) polygon which intersects the 
edges orthogonally at their centers. 

Remark 10 We specifically did not use the word Riemannian because it is not 
entirely clear what Riemannian should mean. Natural guesses would be that 
dij > for all directed edges or that all dual volumes are positive. However, 
we chose not to make such a distinction in this paper. 

2.2 Curvature 

In this section we define curvatures of piecewise flat metrized manifolds, many 
of which are the same as those for piecewise flat manifolds described by Regge 
[55] and Cheeger-Miiller-Schrader [S]. Generally, curvature on a piecewise flat 
manifold of dimension n is considered to be concentrated on codimension 2 
simplices, and the curvature at a is equal to the dihedral angle deficit from 2tt 
multiplied by the volume of ct, possibly with a normalization. Cheeger-Miiller- 
Schrader show that, under appropriate convergence of the triangulations, 
such a curvature converges in measure to scalar curvature measure RdV. (In 
fact, Cheeger-Miiller-Schrader prove a much more general result for all Lipschitz- 
Killing curvatures, but we will only consider scalar curvature.) We first define 
curvature for piecewise flat manifolds in dimension 2, which is concentrated at 
vertices. 

Definition 11 Let{M,T,tj be a two-dimensional piecewise flat manifold. Then 
the curvature Ki at a vertex i is equal to 

Ki = 27r - ^7ijfc, 
where ji are the interior angles of the triangles at vertex i. 
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Angles can be calculated from edge lengths using the law of cosines. Note 
that in two dimensions, curvature satisfies a discrete Gauss-Bonnet equation, 



i 

where x is the Euler characteristic. 

In dimension 3, the curvature is concentrated at edges. 

Definition 12 Let {M, T, £) be a three-dimensional piecewise flat manifold. 
Then the edge curvature Kij is 

Kij = 27r — '^^Pij^M iij- 
\ k.t ) 

The dihedral angles can be computed as a function of edge lengths using the 
Euclidean cosine law to get the face angles, and then using the spherical cosine 
law to related the face angles to a dihedral angle. 

There is an interpretation of Kij/£ij in terms of deficits of parallel transla- 
tions around the "bone" (See [3S] for details.) For this reason, one may 
think of Kij/iij as some sort of analogue of sectional curvature or curvature 
operator (see [7]). 

The fact that curvature is concentrated at edges often makes it difficult to 
compare curvatures with functions, which are naturally defined at vertices. For 
this reason, we will try move these curvatures to curvature functions based at 
vertices. 

In the smooth case, the scalar curvature has interesting variation formulas. 
For instance, we may consider the Einstein-Hilbert functional, 

£n{M,g)^ [ RgdVg, 

J M 

where Rg is the scalar curvature and dVg is the Ricmannian volume measure. 
Note that if rt = 2, then the Gauss-Bonnet theorem says that £'H{M,g) = 
2Trx (M) , but otherwise this functional is an interesting one geometrically. A 
well-known calculation (see, for instance, [3 ) shows that if we consider variations 
of the Riemannian metric Sg — h on M", then 

6£H{M,g)[h]^ f giE,h)dV, (2) 

where E = Rij — \Rgij is the Einstein tensor. It follows that critical points of 
this functional satisfy 

R^J - \R9^, = 0. (3) 

Taking the trace of this equation with respect to the metric, we see that, if 
n ^ 2, this implies that 

R^3 = 0, (4) 
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which is the Einstein or Ricci-flat equation. It also makes sense to consider 
either the constrained problem where volume is equal to one, or to consider the 
normalized functional 

gH(M",g) 

where V is the volume. In both cases we find that critical points under a 
conformal variation correspond to metrics satisfying 

R,j = \g^j (5) 

for a constant A. Taking the trace and integrating, we see that 

l gH(M",g) 

We now consider Regge's analogue to the Einstein-Hilbert functional on 
three-dimensional piece wise flat manifolds. 

Definition 13 If (^M^,T,tj is a three-dimensional piecewise flat manifold, the 
Einstein-Hilbert-Regge functional EHTZ is 

id 

where the sum is over all edges {i,j} G E (T) . 

The analogue of the first variation formula ^ is 

^SrnZ (M, T,i)^2n-Y, P^,M■ (8) 

This was proven by Regge [35^ and follows immediately from the Schlafli formula 
(see [30]). By analogy with the smooth case, we define the following. 

Definition 14 A piecewise flat manifold (^M^,T,tj is Ricci flat if 

K,, = 

for all edges {i, j} . It is Einstein with Einstein constant X (zR if 

dV 

— , (9) 



for all edges {i,j} , whe 



is the total volume. 



i,j,k,£ 
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The term on the left of ^ can be made more explicit. Note that 



3V{M,T,£) = 



da 
d_ 
da 



a^V (M, T, () 



V (A/, T, a I) 



a=l 



dV 



so 



£mZ{M, T,£) 



3V{M,T,e) ' 

analogous to the smooth formula ([6|. Furthermore, we can explicitly compute 
for any tetrahedron {i,j, k,£} that 



ijk£ 
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(10) 



For brevity, we omit the proof of ( 10 1 since we will not use it. However, it can be 
proven by a direct computation of the derivatives of volume and of the dihedral 
angle. 

As in the smooth case, studying the Einstein equation is quite difficult. 
Progress can be made by considering only certain variations of the metric. If 
one takes 6g — fg for a function /, we have a conformal variation. Under 
conformal variations, the scalar curvature satisfies 



Since, under this variation, SdV 
variation is 



SR[fg] = {l-n)Af-Rf. (11) 
^fdV, the variation of EH under a conformal 



SSn{M,g) [fg] = 



J [(l-n)A/+(^-l) Rf 



dV 



In particular, (^ — l) i? is the gradient of £H with respect to the L2 (M, dV) 
inner product. We see that critical points of the functional under conformal 
variations correspond to when the scalar curvature is zero. Note that if we 
either (a) restrict to metrics with volume 1 or (b) normalize the functional, 
then we get constant scalar curvature metrics as critical points. The second 
variation of ETC can be calculated from ( 11 ) to be 

S^Sn {M,g) [fg.Jg]-^ 



1 



- 1 



M 



(l-n)/A/+(^-l) Rf 

("-i)|v/r + Q-i)i?/^ 



dV 
dV. 
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The second variation can be used to check to see if critical points are rigid, i.e., 
if there is a family of deformations of critical metrics. 

The discrete formulation is motivated by the work of Cooper and Rivin [12] , 
who looked at the sphere packing case. The goal is to formulate a conformal 
theory in the piecewise flat setting which allows simple variation formulas as in 
the smooth setting. First, we define the scalar curvature. 

Definition 15 The scalar curvature K of a three-dimensional piecewise flat, 
metrized manifold (^M^,T,d) is the function on the vertices defined by 



-ft'i = ^ I 27r - ^ f3ij^ke | dij. 

k,i 

This definition is much more general than the one in [12j, but restricts to 
almost the same definition in the case of sphere packing (see Section [s] for the 
details). This curvature is in many ways analogous to the scalar curvature 
measure RdV on a Riemannian manifold. Note that, unlike the edge curvatures 
Kij^ this curvature depends on the metric, not only the piecewise flat manifold. 
We also note the following important fact. 

Proposition 16 // {AP,T,d) is a three-dimensional piecewise flat, metrized 
manifold, then the Einstein- Hilbert- Reg ge functional can be written 



Proof. Simply do the sum and recall that dij + dji = iij . ■ 

Now let us define conformal structure. The motivation for the definition will 
be seen in Theorems [23] and [24] and we will see some examples in Section [5] 
The reader may want to recall Definition [6] 

Definition 17 A conformal structure C {M,T,U) on a triangulated manifold 
{M,T) on an open set U C V (T)* is a smooth map 

C{M,T,U) : U ^mzt{M,T) 

such that if d — C {M,T, U) [/] then for each G -B+ and k CzV, 

^f^ 
and 

— ^ = 
dfk 

if k ^ i and k ^ j . 

Notation 18 Often we will suppress the U and simply refer to the domain of 
the conformal structure C {M, T) . 



- d - 
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We can also define a conformal variation. 



Definition 19 A conformal variation of a piecewise flat, metrized manifold 
(Af , T, J) is a smooth curve f : {—e,e) — > V {T)* such that there exists a con- 
formal structure C {M, T, U) with f (— £, e) d U and f (0) = d. We call such a 
conformal structure an extension of the conformal variation. 

An important point is that if we have a conformal structure or conformal 
variation, quantities such as -jrr- make sense. We will usually try to make 

statements in terms of ^ in order to reveal the appearance of discrete Lapla- 
cians, however sometimes it will be more convenient to express terms as partial 
derivatives. We note that a conformal variation is essentially independent of 
the extension in the following sense. 

Proposition 20 Under a conformal variation f (t) of (M, T, (J), we have at 
t = that 

dt '' ~ dt ^ dt ■ 

In particular, at t = 0, for a given ^ (0) , the variation of the length is indepen- 
dent of the extension. 

Proof. From the definition of conformal structure, we have 

_ dijj dfi dljj dfj 

dt ~ dfi dt dfj dt 

dfi dfj 
^d,j— + d,, — . 



Notation 21 In the sequel, when we suppose a conformal variation, it will he 

dt 



understood that quantities such as ^ are evaluated at t — 0, though not stated. 



There are often more than one extension to a conformal variation. For 
instance, for a triangle {1, 2, 3} , we may extend the metric defined by — ^ 
for all (z, j) S E+ to several families where fi (t) = txi, such as 

dij (t) = ^ exp (tXi) , 



which corresponds to a circle packing conformal structure (see Section 5.1 1, and 



dij (t) = ^exp ^ {x, + Xj) 



which corresponds to a perpendicular bisector conformal structure (see Section 
53|. 
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Remark 22 Often a conformal structure will be generated from a base metric, 
much the same way a conformal class on a Riemannian manifolds can be de- 
scribed as the equivalence class of metrics e-^ go, where f is a function on the 
manifold and go is the base Riemannian metric. However, we have not defined 
it thus and, in general, one must be careful how the structures are defined if 
one wishes to partition all piecewise fiat manifolds into conformal classes. We 
do not attempt this here, though there is a straightforward way to do this for 
perpendicular bisector conformal structures seen in Section\5.3\ 



In two dimensions, the fact that curvatures arise from conformal variations 
of a functional is not obvious, but can be proven. 

Theorem 23 Fix a conformal structure C (^M^ ,T,U^ on a two dimensional 
triangulated manifold and suppose that U is simply connected. Then there is a 
functional F : U M. such that 



dF 



for each i £ V (T) . Furthermore, the second variation of the functional under a 
conformal variation f (t) can be expressed as 



(fF_ 
df^ 



9 /-^ 



dt 



df^ 
dt 



dt^ 



(12) 



This sort of formulation of the prescribed curvature problem in a variational 
framework has been studied by many people. See, for instance, [T3] [33] [TT] [5] 
IS] [H]- However, no source to date has unified the theorem in the way of 
Theorem m 

In three dimensions, the Einstein-Hilbert-Regge functional is a natural one 
to consider. 



Theorem 24 For any conformal variation f (t) of a three dimensional, piece- 
wise flat, metrized manifold [M, T, cF) , we have. 



di 

d^ 
dt^ 



snn{M,T,t{fm = Y,K^. 

i 

£nn{M,Tj{f{t)))^Y.T. 



dt 



(13) 



0.. o/- 



df^_dfA 

dt dt J ^ ' 



dt^ 



where 



9y 



ddij ddji 
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Thus a critical point of STiTZ corresponds to when Ki = for all i and at a 
critical metric, 

^£KR(M,r,.(/W))^i:i:(|-|^A-,) (f -f)^ (15) 



Furthermore, at a critical point for the general ETCR- {M, T, £) , we must have 
Kij = (see W^), and hence here we have 



Theorem 24 motivates the definition of constant scalar curvature metrics, as 
seen from the following. 

Corollary 25 For any conformal structure of (Af , T) , we have 

-§-Enn{M,T,i{f))^K,, 

oft 

^V{M,T,i{f)) = V., 
where ^ 



The proof will be given in Section |4.2[ Corollary [25] motivates the following 
definition. 

Definition 26 A three-dim,ensional piecewise flat, metrized manifold (^M^,T,d) 
is has constant scalar curvature A if 

= AFi 

for all vertices i. 

It is not hard to see that 



^ F. = 3V. 



Summing both sides of the constant scalar curvature equation, we see that 

1 EWl (Af, T, i) 



A = 

Note that 



3 V(A//,T,^) 



T/. _ ^ _ \- 

'~ dfr dhp'' 
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and so on an Einstein manifold, which would satisfy 



for each edge, we have that 

3 3 

We have just proved the following. 

Proposition 27 // {AP,T,d) is a three-dimensional piecewise flat, metrized 
manifold which is Einstein, then it has constant scalar curvature. 



There is a second variation formula for conformal variations of ETLTZ/V 



1/3 

at Einstein manifolds, but for brevity we omit it since it requires the calculation 
of Iji. With the results from this paper, it is straightforward to calculate these 
derivatives. 



3 Variations of angles 

In the rest of this paper, we will use 5 to denote the differential. 
3.1 Two dimensions 

In this section we will compute the derivative of an angle under a certain 
variation of lengths. Consider the Euclidean triangle determined by lengths 
(^12, ^13, ^23) with vertices {vi^v^^vz} and also the triangle determined by lengths 
(^12: ^13 + '5^13, -^23 + ^^23), say with vertices {wi, W2, ^^3}, under the important 
assumptions: 

5h2 = 0, (17) 

543 = ^31^/3, (18) 
5^3 = ^32^/3, (19) 

where d is a metric on {1, 2, 3} inducing lengths £. 

Draw the arc representing £14 (57i_24, which goes through vertex vj, and 
intersects the segment viv'^. Call this edge E. It has endpoints and w. See 
Figures [1] and [2] for cases when the center C123 is inside W1W2W3 and when it is 
outside W1U2W3. 

Proposition 28 The points C123, 1^3, and v'^ lie on a line. I.e., 5v^ is parallel 

to Vz - Ci23- 
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Proof. Notice that 

^ (^?3) = S [iv3 - V,) ■ - Vi)] = 2iv3- Vi) ■ SV3 

but also 

<5 (4) =2^13^31(5/3, 



SO 



Similarly, 



OV3 ■ — = dsidfs. 

«13 



. V3-V2 „ 
OV3 ■ — -„ = 0.320 j^. 



t23 

Similarly, the vector 113 — C123 satisfies 



{v3 - C123) ■ = dsi, 

*13 

. . V3-V2 

{V3 - C123) • -„ = 032 

t23 
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Vl V2 

Figure 2: Variation of a triangle where C123 is outside the triangle. 



and so we see that Sv^ = {v^ — C123) Sf^. m 

This is essentially the same proof given by Thurston [44] and Marden- Rodin 

m- 

Consider the triangle V3WV31. This is a right triangle with right angle at w 
since viw is a radius of the circle containing E. Since the angle of E with V1V3 is 



also a right angle, together with Proposition 28 it follows that V3WV3' is similar 



to the right triangle C123C13W3. Using the similar triangles, we get that 

^13 ('^71,23) _ ^13,2 ,„„s 

«13 " d3i ■ ^^"^ 

So 

t^7l,23 __ fel3,2 

5/3 " hs ■ 
This leads to the following theorem. 

Theorem 29 For variations of the lengths of a Euclidean triangle of the type 
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77pP (where Sfs is arbitrary), we have 



djl.23 


hl3,2 




h3 


9^2,13 


^■23,1 


dh 


^23 


973,12 


/ll3,2 ^23,1 


dh 


^13 ^23 



Proof. We have already proven the first two equahties. The last follows from 
the fact that in a Euclidean triangle, 71,23 + 72,13 + 73,12 = tt. ■ 

Remark 30 Special cases of this theorem appear in f23 \/\17}jf l 8 f f 1 9 f and less 
refined versions (where only signs and not explicit values are computed for the 
derivatives) appear in many other places, including \44^i39fl29fl4 

3.2 Three dimensions 

Now consider a tetrahedron {1, 2, 3, 4}. Similarly, we will need variations of the 
form 

S£i2 = 0, (21) 

S£i3 = 0, (22) 

S£23 = 0, (23) 

5^14 = dii6f4, (24) 

«24 = d42^/4, (25) 
«34 = ^43^/4, (26) 

where d is a metric on {1, 2, 3, 4}. For convenience, we embed the vertices of the 
tetrahedron as {vi,V2,V3,V4} so that it has the correct edge lengths and such 
that vi is at the origin, V2 is along the positive cc-axis, ^3 is in the xy-plane, 
and V4 is above the xy-plane. We will need some additional points. We let v'4 
be the new vertex 4 gotten by taking lengths £ij + S£ij (remembering that, for 
instance, d£i2 = 0) and embedding this tetrahedron as viV2V3v'^, with v'4 above 
the a;y-plane. We also need v'^ 12 which is the point on the plane ViV2v'^ which 
makes viV2v'^i2 ^ triangle congruent to V1V2V4. Also, we have the point v'^^, 
which is the point on the line W1W4 which is a distance £14 from vi. See Figure 

El 

We first observe that the right tetrahedron 12V4 iv'^ is similar to the 
tetrahedron C1234C124C14W4. This is because C1234, V4, and v'^ are colinear (the 



proof is exactly analogous to the proof of Proposition 28 and is thus omitted). 
This implies that 

(5^14)^ _ 5 (^14<57l,24) (^14Sin7i.24 ((5/3l2,34)) 

9/114,2^124.3 
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V2 

Figure 3: Variation of a tetrahedron. 

We conclude that 

^124,3^14, 2 _ (^71,24 ^ / '^/3l2,34 \ 

sin 71^24^14 f^^w / 

dil V '5^14 / 

using Theorem [29] We thus get 

/ll24,3 _ <^/3l2,34 .27^ 

d41 sin 71^24^14 <5^14 

Furthermore, if ai is the solid angle at vertex vi, then we get that ll^^Sai is 
approximately the sum of the areas of two triangles on the sphere centered at vi 
of radius £14. One triangle has vertices W4 1, W4, and a point on the x-axis which 
we will call b. See Figure g] This tr iangle can be divided into two spherical 
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V2 % 



Figure 4: Angle variation setup for a tetrahedron. 



triangles, U4 ^v'^ 12V4 and v'^ i2V4:b. Each of these triangles has a right angle at 
i'4 Note that the first triangle has area which vanishes to higher order, so 
up to first order, the area is the area of the second triangle. This triangle has 
a right angle at v'^ angle 5/3i2,34 at 6, and some other angle at V4, say 

^ — 7, where 7 is small. It is easy to see that the side v'^ in the spherical 
triangle has length £146^1,24 and the side W4 i2'^4 has length ^14 sin 7i.24(5/3i2,34 
to first order (the error is like (^/3i2,34)^)- We may compute 7 since we know by 



spherical trigonometry that 

cos (I -7) 



tan ((5/3i2,34 sin 7i_24) 



tan 71^24 

We look at this asymptotically where 5/3i2,34 is small, and see that 

<5/3l2,34 sin 71,24 .r, 
7 = 7 = 0P12,34 cos 71,24 

tan7i 24 
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and thus the area of the triangle is 

('^/3l2,34 + f + f ~ '^/?12,34 COS 71,24 " Tt) 
= f u^Pl2,3i (1 - COS 71,24) • 

Hence we have that 

5ai = 5/3i2,34 + ^/3l3,24 + <5A4,23 

impUes that 

(5/312,34 (1 - COS 71,24) + ^/3l3,24 (1 - COS 71,34) = (5/3i2,34 + ^/3l3,24 + ^/3l4,23 



or 



'5/3l4,23 = -COS7i,24<5/3l2,34 " COS 7i,345/3i3,24 

^124,3 c-/) 

l-j -. J—OI14- COS. . 

d4iSin7i,24«14 041 sin 71,34^14 



^124,3 c-/) /ll34,2 
= -COS7i,24;^ -. ■ ^^«14 - COS7i,34-5 : ' ^0«14 



d. 



(cot 71,24/1124,3 + cot 71,34/1134,2) d£i4 



41*14 



1 ^ u I. , u u C^12/ll24,3 C?13/ll34,2 1 c-„ 

'Il4,2'il24,3 + 'Il4,3ftl34,2 ^ : 0<;i4- 



dudiilu V ' ' ' ' sin 71,24 sin 71,34 

Recalhng that 



il24.3 <5/3l2,34 



we get 



That is, 



^41 sin 71,24^14 %4 



d41 ((ii45/3i4,23 + (il2'5/3l2,34 + <^13<5/3l3,24) 

= -7 — (/ll4,2/il24,3 + /ll4, 3/1134,2) ^^14 
«14 

= 2^5^14 = 2^d4i<5/4. (28) 

«14 «14 



<^14<5/3l4,23 + t^l2<5/3l2,34 + rfl3'5/3l3,24 = 2 ^]'^''^^ 6/4. 

l\4 



Furthermore, the Schlafli formula implies that in the tetrahedron, 

42^/312 + 43^/3l3 + 44^/3l4 + ^23<5/323 + ^24^/324 + 44%4 = 0, 

and so 

= ((ii4(5/3i4 + di26f3i2 + disSPis) + {d2id(3i2 + rf23<5/323 + (^241^/324) 

+ (^31^/313 + £^32^/323 + (^34^/334) + (c!4l5/3l4 + ^42^4 + c!43<5/334) 
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The first three terms are all of the form ( 28 1 and the last is different (because 
the lengths of edges around vertex 4 are changing). Thus, 

(341*14 "42«24 «43«34 



A 



14,23 



A 



24,13 



A: 



t-24 



34,12 



'^34 



5h. 



We have just proven the following theorem. 



Theorem 31 Under variations of the form (21)-(26), we have 



(^24(5/324 
^34(^/334 

diiS(3i4 



' 0^23 '5/523 
■ (^32 (5/323 
- C?42 (5/324 



■ d2lSPi2 

dsiSPia 
disSPsi 



„ ^24,13 cr 

2— S,U, 

t24 

„^34,12 cr 

2— (5/4, 

t34 

„^14,23 



r) ^24,13 _ ^34,12 
^4 44 



5/4. 



We actually derived a finer result, with explicit computation of the variations 
of individual dihedral angles. However, the result in this form is more compactly 
stated and all we will use in the remainder of this paper. Also, this result 
was derived in |17| for the specific case of sphere-packing configurations of a 
tetrahedron. Less precise results along these lines were examined in the sphere- 
packing case in [12] [37] as well. We will go into more detail later. 



4 Curvature variations 



4.1 Two dimensions 



Discrete curvature in two dimensions has been well studied, with curvature as 
in Definition [11] Theorem [29] has the following implication. 



Theorem 32 Let f (<) be a conformal variation of (M^, T, d) Then for t = 0, 

(29) 



dK, 
dt 



f* 



df^^df^ 
dt dt 
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Proof. Wc compute 



dt 




which imphes the result. 



Remark 33 We will express many of the variation results as in Theorem 32 
instead of in terms of in order to emphasize the presence of the Laplacian. 
Notice that the formula (29) has the form 



dJU 
dt 



A 



for an appropriate definition of the Laplacian A. We will comment more on 
this in Section\^ 



The formulas from Theorem 
giving the proof of Theorem ^ 
Proof of Theorem 1231 F will be defined as 



29 also imply that the curvatures are variational 



F{f) 



f 



fo 



for the 1-form 



UJ ■ 



To ensure that the integral is independent of path, we need that ui is closed. 
Since we are in a simply connected domain, we need only check that 



cm 



a/. 



for i 7^ j. We can compute these derivative explicitly, and they are 



if i and j share an edge and zero otherwise. ■ 

This formulation of the prescribed curvature problem in a variational frame- 
work has been studied by many people. See, for instance, [13] [33] [TT] [S] [H] 
|21j . many of which derive the functional in precisely the same way. Certain 
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conformal variations in the discrete setting have been proposed by Rocek and 
Williams in the context of Regge calculus [3S] , Thurston in the setting of circle 
patterns [44] (see also other work on circle packing, e.g., |l3]), and Luo [28] 
(see also [H]). Our current setting takes each of these definitions and proofs 
as special cases (see Section [s]) . In many of these papers it is shown that the 
largest "reasonable" domain for the f's is simply connected. The advantage of 



Theorem 23 is that it works for a more general class of conformal variations. In 
addition, we have a geometric descr 
from most of these previous works. 



addition, we have a geometric description of the derivatives which is absent 



4.2 Three dimensions and the Einstein-Hilbert-Regge func- 
tional 

We could follow a similar method to that in Theorem [23] to prove that there are 
three-dimensional curvatures which are variational. However, we will present 
this fact in a different way by using the Einstein-Hilbert-Regge functional from 
Definition [131 

Recall Definition [15] for scalar curvatures in dimension three. The definition 
is first motivated by seeing how the Schlafii formula decomposes as sums around 
vertices. The Schlafii formula on a tetrahedron is 



=0, 



where the sum is over all edges in the tetrahedron. It can be written as 

= (di2 + d2l) Spi2 + {dl3 + d3i) SPi3 + {di4 + d4l) S(3i4 
+ (^23 + d-i2) SP23 + (^24 + d42) 5(324 + (^34 + ^43) 1^/334 
= (c!l2<5/3l2 + (il3<5/3l3 + di45l3i4) + (rf21<5/3l2 + d23SP23 + ^24-5/324) 
+ {d3iS/3i3 + d32<5/323 + ^34(^/334) + id4lSPi4 + (i42%4 + rf43'5/334) , 

giving the vertex breakdown motivating the curvature formula. The Schlafii 
formula allows an easy computation of first derivatives of the Einstein-Hilbert- 
Regge functional on a triangulation of a closed manifold, giving 



k,e 



k,e 



3 
K, 



(30) 



To compute the second variation, we need the variation of Ki. Using Theo- 



rem 31 we can express the derivatives of curvature. 
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Theorem 34 Let f (t) be a conformal vari ation of {M^,T,d) . Th en, 



dt 



dt dt 



J2 I 27r-I]^^J'fe^ I Jt'^'i 



(31) 



E(2 



dfj dfi 
dt dt 



dh 
dt 



wht 



ddij ddji 



We note that the first term in (31 ) is the Laplacian operator studied by the 
author in [T5] . 
Proof. We compute 



dt ^ Hf "-^ ^ 



P- ■ 



dfj dfi 
dt dt 



( Pi3M~ndij + f^ikje-rrdik + Pu,jk-T,du 



dt 

ij ( dfj dfi 
dt ~ ~di 



^ \ k,e J 



dt 



Furthermore, if there is a conformal structure, then 



dfj ' dfj df, " df, dfj " df. 



d diij d diij ddji 



Thus 



ddij ddji d'^ii 



dfj dfi df^dfj 
is symmetric, i.e., qij = qji. This also implies that 

ddij diij ddji 
dft df, dfi 



Thus 



d , \dfi ^ dfj 



The result follows. ■ 

We can compute the variations of the Einstein-Hilbert-Regge functional. 
Proof of Theorem 24 , Using ( 30 ) , we see immediately that 

df^ 



d 
dt 



£mZ{M,T,£if{t)))^J2^^ 



dt 
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Using Theorem 34 we compute that 



dt dt 



EE 



df^ 



dfj _ dfi\ dii^ 
dt dt I dt 



dt 



dt^ 



£* 



dfj df, 
dt dt 



E^^ 



dt 



dt^ 



The rest of the theorem follows immediately from (30) and Regge's variation 
formula Q. ■ 

In order to get better control of qij , we will look at special conformal struc- 
tures in Section [S] 

Finally we can complete the proof of Corollary [25] 
Proof of Corollary |25[ We see from Figure [3] that we must have that if we 
fix /i, /2, /a and let vary, then 

^Vi234 = ^^124^14 sin 71, 24^/3l2,34 + ^^134^14 siu 71, 34(5/3i3, 24 



1 



7^23444 sin 73,24<5/323, 



{Ai2ih 



124,3 



^134^134,2 



234 "■2344^ 



by ( 27 1 , which implies that 



where 



— V{T,£{f)) = V, 



ijk 



3,k,e 



and the sum is over all tetrahedra containing i and all faces in those tetrahedra 
containing i. ■ 



5 Examples of conformal structures 

In this section we place previously studied geometric structures into the frame- 
work of conformal structures. 
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5.1 Circle and sphere packing 



The case of circle packing and sphere packing is when edge lengths arise from 
spheres centered at the vertices which are externally tangent to each other. In 
this case, there are positive weights corresponding to the radii and £ij — Vi+rj. 
In two dimensions, circle packings have been considered in a number of contexts; 
see Stephenson's monograph [43] for an overview. In three dimensions, this case 
was considered by Cooper- Rivin |12j . They noticed, in particular, that for a 
sphere packing, one can rewrite the Schlafli formula in the following way: 



= (ri + ra) (5/3i2 + {ri + rg) S(3i3 + {n + ri) SPu 
+ {r2 + rs) 6/323 + {^2 + ta) (5/324 + {rz + ri) Sp^i 

= r-i ((5/3i2 + <5/3i3 + 6l3u) + r-2 (<5/3i2 + <5/323 + '5/324) 
+ (^/3i3 + (5/323 + ^(3^^) + {6(3u + 6(324 + SPsa) 

= ri6ai + r26a2 + r^Sas + r46a4, 

where ai is the solid angle at vertex i, and thus 6ai — 6 (/5i2 + P13 + /3i4 — tt) = 
<5/3i2 + (5/3i3 + (5/3i4. They used this to motivate the definition of scalar curvature 
as 47r — ^ ai where the sum is over all tetrahedra containing i as a vertex. From 
our setting, we would define the scalar curvature measure instead as 

V J ] \ k,i 

where in the right side, the first sum is over all edges incident on i and the 
second sum is the sum over all tetrahedra containing as an edge. The 

second equality can be easily derived using the Euler characteristic and area 
formula of the sphere centered at vertex i. 

To match this to our setting, we see that we must take fi = logr^ and 
dij = Ti, since 

^ij — r^ Tj 

-?r^n = n. 

Thus we have the following conformal structure. 

Definition 35 The circle/sphere packing conformal structure, {M,T) , is 
the map defining 

dij = e-^' 

for every oriented edge in (T) restricted to an appropriate domain of f G 
V* (T) . 

In two dimensions, the triangle inequality is automatically satisfied, and 
so the domain is all of V* (T) . In three dimensions, there is an additional 
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condition that the square volumes of three-dimensional simplices (as defined by 
the Cayley-Menger determinant formula) are positive. This is discussed in some 
detail in [17| [E]. 

We see that the formulas in Section [3?2| correspond to 

dai „ ^14,23 
^— rir4 = 2— . 

This is the same formula derived by the author in fl7j . 



5.2 Fixed intersection angles/inversive distance 

There is a more general case of circles or spheres with fixed intersection angles, 
originally considered by Thurston |44j . Here we parametrize lengths by two 
parameters, radii rj and inversive distances 77^. The inversive distance (see, 
for instance, [H]) is like the cosine of the supplement of the intersection angle, 
defined so that 

il^r^ + r^ + 2nr,7^,,. 

We will use this formula to parametrize the lengths by the radii with inversive 
distances fixed. It essentially corresponds to having circles at the vertices of 
radius and intersecting at angle arccos (— r/ij) . If tjij is not between —1 and 
1, then the circles may not intersect, but this is not a problem for the theory. 
There is always a circle orthogonal to these circles, and we take the center of the 
triangle to be the center of this orthocircle. (Note, it is possible that this circle 
does not have real radius, but the center is still well defined using the algebra 
of circles given in [32,.) We then find that 

d^,^':i^^:i±^. (32) 

For a path in the r variables, we compute 

„ „ . . dri drj 

dr'~^''n dt ^ dt ■ 

Thus we see that fi = logr^, giving the fixed inversive distance conformal class. 

Definition 36 For a given rj G E [T)* , the fixed inversive distance conformal 
structure, C^^ (A/, T, 77) , is the conformal structure described by the map 



% (/) 



where 

is the length, when restricted to a proper domain in V (T)* 
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Note that there are some restrictions on the domain which may be quite 
complicated, including the triangle inequality. However, it has been found that 
in two dimensions, if rjij > for all ij G E (T)* , then the domain is simply 
connected. This was initially shown for < rjij < 1 by Thurston ([33] [53]) and 
the additional cases were proven recently by Guo |21j . 

We see that 

ddi-j ddi 



dfj dr 

.2 
j 



£3. 



We finally get 

dK, ^ 
dt 

and the second variation of the Einstein-Hilbert-Regge functional is 




i j^i 




dfj dfi 
dt dt 



Note that in the case that r/ij — 1 , corresponding to sphere packing, the second 
term is zero. In general, for spheres with intersection we have rjij < 1 and for 
spheres which do not intersect we have r]ij > 1 and so in each case the term 
with edge curvatures has a particular sign. 

5.3 Perpendicular bisectors 

Here we give the conformal structure proposed by Rocek- Williams [38], Luo 
[28] , and Pinkall-Schroeder-Springborn [42]. This structure has also been found 
in the numerical analysis literature on approximations of the Laplacian in the 
context of the box method (see, e.g., [5^ and [33])- Take 



where Lij are fixed lengths. We see that, given a path in the space of u variables. 



If we take 



/. . — /. . 

dt '' ~ 



dui 
~dt 

I. 



2 

2ui 



duj 
~dt 
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^ dt 



then 

dt '' ~ dt 

We notice that the duals to the edges intersect the edges at their midpoints, 
which is why we call this the perpendicular bisector conformal structure fol- 
lowing [55]. It can be proven inductively that the center of any simplex is the 

center of the sphere circumscribing that simplex. 



Definition 37 Let L Cz E (T) be such that {M, T, L) is a piecewise flat man- 
ifold. The perpendicular bisector conformal structure, C^^ (M, T, L) , is the 
conformal structure determined by 



y 2 ■ 



when restricted to an appropriate domain. 



Since 



(m, T, d 



is a piecewise flat, metrized manifold, this conformal 
structure exists for fi close to 0. However, the largest possible domain must 
satisfy a number of inequalities. 
We see that 

_ dd^ _ 1 
'^'^ " dfj ~ A ''- 

Thus, in three dimensions the variation of curvature is 



dt 



E 



1 



dt 



dt 



dh 
dt 



We get that 



dt 



^snn (T, ^ (/ m = E E ( # - s 



E^^ 



dh 

dt 



df^ 
dt 





d^f^' 




dt^ 



6 The discrete Laplacian and the second varia- 
tion 

6.1 Laplacians 

The relationship between the second variation of the functionals presented here 
and the Laplacian is the main reason we describe these variations as conformal. 
The standard Laplacian is defined as follows. 
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Definition 38 Let (M, T, d) be a piecewise flat, metrized manifold. The dis- 
crete Laplacian A is an operator V (T)* — > V (T)* defined by 



for each vertex i, where £*j is the dual length defined appropriately (see Defini- 
tions^ and^ and I19f for the general case). 

These can be considered Laplacians on the graph of the 1-skeleton with 
edges weighted by j^. (For more on Laplacians on graphs, see [TO]-) This is a 
very natural choice of Laplacian, arising, for instance, by considering another 
function tp on the vertices, and defining the Laplacian weakly as 

I I, J 

for all choices of iJj, where Vij is the volume associated to an edge, defined by 

V- ■ — -/*■/■ ■ 
n 

where n is the dimension. This is an analogue of the definition of the smooth 
Laplacian on a closed manifold as the operator such that 

A(t)^p dV = - J dV 

for all smooth functions ip. Another interesting observation about the Laplacian 

I* ■ . 

is that the weights are very much like conductances, in that they are inversely 
proportional to length and directly proportional to cross-sectional area if one 
considers current through wires located at the edges of the triangulation. 

Laplacians of this geometric form have been studied for some time. The most 
well-known is the "cotan formula" for a Laplacian on a planar triangulation. If 



one considers the perpendicular bisector formulation of Section |5.3| on a planar 
domain or surface, one finds that 1*^ = iij {cot jk,ij + cot ji^ij) . It turns out 
that this is precisely the finite element approximation of the Laplacian, as first 
computed by Duffin [TS]. The cotan formula has been well-studied both in 
regards to approximation of the Laplacian on domains and approximation of 
the Laplacian on surfaces for computing minimal surfaces and bending energies. 
See, e.g., [SS] [23] [IS] [S]. In addition, Laplacians have appeared in the 
study of circle packings. In fact, to our knowledge, the first observation that 
variations of angles are related to dual lengths dates to Z. He [23] in the circle 
packing setting, where it was used for constructing a Laplacian. Further work 
in two dimensions in the setting of circle packings and circle diagrams with fixed 
inversive distance which connects angle variations with Laplacians can be found 
in [15] [9] [20] [21] . An interesting study of possible Laplacians from a axiomatic 
development can be found in [47]. 
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6.2 Properties of the Laplacian 

There are two properties of the smooth Laplacian which are desirable to have 
in a discrete Laplacian: 

1. A is a negative seniidefinite operator with zero eigenspace corresponding 
exactly to constant functions (0 is a constant function if there exists c G M 
such that = € for all i £ V). 

2. A satisfies the weak maximum principle, i.e., for any (f> £ V (T)* , if 0^ = 
mini 0i and 0m — max^ then A(/i,„ > and l\(j)M < 0. 

Note that the definition of the Laplacian ensures that the constant functions 
are in the nuUspace. The second property is implied by (*j > for all edges 
{i,j} ■ Furthermore, we shall show that the strict inequality £*j > implies the 
property [l] The Laplacian is a symmetric operator, and so it has a full set of 
eigenvalues. If A is an eigenvalue with eigenvector 0, then 

= E 7^ - ■ 

3 ■' 

We see that 




and so we see immediately that if i*^ > 0, then A < 0. Furthermore, if the 
inequality is strict, then A implies that 0^ = for every edge. On a 
connected manifold, this implies that (j) is constant. This type of Laplacian 
has good numerical properties and for this reason numerical analysts are often 
interested in using such a Laplacian for numerical approximation of PDE. (For 
instance, see [25].) 

In two dimensions, the property 1*^ > is a weighted Delaunay condition 
[TO] . Note that the argument in the previous paragraph shows that this condi- 
tion implies that the Laplacian is negative semidefinite, but it may have a larger 
nuUspace than just the constant functions. Often in triangulations of the plane, 
one gets around the fact that the inequality is not strict by removing edges 
with the property that £*j — and replacing the triangulation with a polygo- 
nalization. In the manifold case, this could potentially introduce curvature to 
the inside of the polygons, so we do not pursue this direction. It is not known 
whether a given piecewise fiat, metrized manifold (M^, T, d) can be transformed 
to another piecewise fiat, metrized manifold [M^ ,T' ,d') which is weighted De- 
launay such that the two induced piecewise fiat manifolds are isometric in a 
reasonable sense. This is true for the perpendicular bisector conformal struc- 
ture, which corresponds to finding Delaunay triangulations (see [36] and [B]). In 
three dimensions, the property i*j > is not equivalent to a weighted Delaunay 
condition, and much less is known about the existence of such metrics. How- 
ever, the geometric description of the Laplacian ensures that if all the centers 
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of the highest dimensional simplices are inside those simplices, then the second 
property is satisfied (some call this property "well-centered," see [H]). 

The first property is certainly weaker. There are a number of instances 
when one can prove the first property without the second property being true. 
For instance, for a metric in a two-dimensional perpendicular bisector confor- 
mal structure, we see that the induced Laplacian is precisely the finite element 
Laplacian. This Laplacian always satisfies the first property, but only satisfies 
the second if it is Delaunay (see [36] for a proof) . We state a proposition sum- 
marizing the known conditions which ensure the first property. The following 
proposition is an amalgam of known results. 

Proposition 39 Let {M^,T,d) be a piecewise flat, metrized manifold. The 
discrete Laplacian is a negative semidefinite operator with zero eigenspace cor- 
responding exactly to constant functions if any of the following are satisfied: 

1. £*j > for all edges {i,j} e E (T) . 

2. n — 2 and the triangulation is in C^^ (Af , T,ri) , a fixed inversive distance 
conformal structure, with rjij > for all {i,j} G E (T) . 

3. n = 2 and d^J > for all {i,j) G E+ (T) . 

4- n — 2 and (M, T, d) is C^^ {M, T, L) , a perpendicular bisector conformal 
structure, for some L. 

5. 71 = 2 and for each triangle isometrically embedded in the plane as ViVjVk, 
the center Cijk is contained within the circumcircle. 

6. n — 3 and (M, T, d) is in {M, T) , the sphere packing conformal struc- 
ture. 

Proof. The proofs follow from a number of results from the literature. The fact 
that ([!]) implies definiteness is well known in the numerical analysis community 
and proven in the discussion before the statement of the proposition. The fact 
that ([2]) implies definiteness was proven for < rjij < 1 by Thurston [ll] and 
Marden-Rodin 29J and the general case of (|2| was proven by Guo [5T]. In fact, 
using ( [32| ), one easily sees that ^ implies (3 1, and the fact that ^ implies 
definiteness is in [TS]. We believe ([3| implies (5]), though we have not verified 
the proof since there is a direct proof for ^ . (|4| implies definiteness was shown 
by Rivin [36 . Also, for ^ , the center is the circumcenter and thus ( 4 1 implies 
([5|. The fact that ([5| implies definiteness is in [5D]. The fact that (5 1 implies 
definiteness follows easily from the definiteness of the related matrix in the 
Appendix from [TH] (also in [T^] and [37j). ■ 

Note that Proposition [39] only covers a small subset of the cases one might 
be interested in. It is of interest that ([2]) -([6]) are all proven by proving the 
definiteness on a single simplex and then extrapolating to the entire complex, 
though ([1]) and takes the global structure into account. In light of ^ and ([s]), 
it may be surprising that the same are not true, in general, for n — 3 (one can 
consider tetrahedra which are close to flat). It would be of interest to know a 
condition similar to ^ which implies definiteness for n ~ 3. 
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6.3 Convexity and rigidity of curvature functionals 

We can use our analysis of the Laplacian to attack two questions about curvature 
functionals: 

Ql. Are the functionals convex? 
Q2. Are critical points rigid? 

The first question is more difficult, but if we take first order variations of fi 
in two dimensions (i.e., = 0), then we have the following theorem. 



Theorem 40 ( [H] |29] [2T] [28] ) The function F described in Theorem^ 
convex on the image of following conformal structures: 

1. C^^ {AP,T,r]), with ri,j > for all {i,j} e E (T) . 

2. C^^ {AP,T,L), for some L. 



Proof. The proof follows immediately from Theorem 23 and Proposition 39 



This theorem was previously proven by combining theorems of the articles listed. 
■ 

In three dimensions, this question is far more complex, much like in the 
smooth case, due to the presence of a reaction term. However, we do have the 
following result. 

Theorem 41 The Einstein-Hilbert-Regge functional is convex on the following 
sets: 

1. Metrics in the image of the conformal structure (^M^,T^ with Ki > 
for allieV (T) . 

2. Metrics in the image of any conformal structure of (^M^,T^ which satisfy 
£*j - IqijKij > for each {i, j} € E (T) and K, > for all i e V (T) . 

We note that (1) is not a special case of (2). In case (1) we have that qij — 
but do not require £*j > 0. We also note that a special case of (2) is a metric in 
the image of C^' {M^,T, rj) with 1 < and K,, > for all {i,j} e E {T) . 
Proof. Recall the variation formula from Theorem 24 As already remarked, 
in case (1) we have qij = 0. Together with case ^ in Proposition 39 the case is 
proven. Case (2) can be proven by essentially the same argument used to prove 
Proposition |39] part ([T|. ■ 

The second question above asks about rigidity, which we can define thus. 

Definition 42 A piecewise flat, metrized manifold (M, T, d) is rigid with re- 
spect to conformal variations if there is no conformal variation f (t) such that 
{M,T, d (f (t))) is fixed other than the trivial variation which scales the edge 
lengths uniformly (in Riemannian geometry, this is called a homothety). 
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Since we have functionals of (M, T, d) in two and three dimensions, we have 



the following immediate consequences of Theorems 23 and [24] together with 
Proposition [39) 

Theorem 43 A two-dimensional piecewise flat, metrized manifold {^M'^,T,d^ 
with curvature zero (i.e., Ki — for all i € V (T) ) is rigid with respect to any 
conformal variations if it satisfies {^-(^ in Proposition \3£\ 

Theorem 44 A three-dimensional piecewise flat, metrized manifold (^M^,T,d) 
which is Ricci flat is rigid with respect to any conformal variations if it satisfies 
|Ip or ^ in Proposition 39 



Note that these statements are analogous to a theorem of Obata in the 
smooth category. 

Acknowledgement 45 This work benefited from discussions with Mauro Car- 
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